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' In this article we show that boundary conditions can be treated as Lagrangian and 

Hamiltonian constraints. Using the Dirac method, we find that boundary conditions 
are equivalent to an infinite chain of second class constraints which is a new feature in 
the context of constrained systems. Constructing the Dirac brackets and the reduced 
phase space structure for different boundary conditions, we show why mode expanding 
and then quantizing a field theory with boundary conditions is the proper way. We 
also show that in a quantized field theory subjected to the mixed boundary conditions, 
the field components are noncommutative. 
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1 Introduction 



It is well-known that to formulate a general classical field theory defined in a box, besides 
the equations of motion one should know the behaviour of the fields on the boundaries, 
boundary conditions. Boundary conditions are usually relations between the fields and their 
various derivatives, including the time derivative, on the boundaries, expected to be held at 
all the times. In Hamiltonian language the boundary conditions are in general functions of 
the fields and their conjugate momenta; hence the field theories subjected to the boundary 
conditions might be understood by the prescription for handling the constrained systems 
proposed by Dirac M. 

In the usual field theory arguments, since boundary conditions are usually linear com- 
binations of fields and their momenta, one can easily impose them on the solutions of the 
equations of motion, and find the final result. But, imposing the boundary conditions in 
some special cases may lead to inconsistencies with the canonical commutation relations 
@, |, |, |, 

In this article, considering the boundary conditions as constraints, we apply the Dirac's 
procedure to this constrained system. Although this idea have been used in || §], the 
problem has new and special features in the context of constrained systems on which, we 
mostly concentrate. 

In the second section, we review the Lagrangian and Hamiltonian constrained systems. 
In section 3, to visualize the seat of boundary conditions we take a toy model and by 
discretizing the model show that boundary conditions are in fact the equations of motion 
for the points at the boundaries so that, when we go to the continuum limit, i.e. , the 
original theory, the acceleration term disappears. In other words boundary conditions are 
Lagrangian constraints which are not consequences of a singular Lagrangian. In section 4, 
going to Hamiltonian picture we study the constraint structure resulting from the boundary 
conditions, and apply it explicitly to some field theories. Implying constraint consistency 
we show that although the Lagrange multiplier is determined, the constraint chain is not 
terminated, This is a new feature in the constrained systems analysis. Exhausting all the 
consistency checks we end up with an infinite constraint chain which all of them are of second 
class, which is another new feature of this constraint structure. Moreover, we construct the 
fundamental Dirac brackets, the Dirac brackets of fields and their conjugate momenta. In 
section 5, by a canonical transformation we go to the Fourier modes in terms of which, the 
constraint chain obtained in the previous section can be easily solved. In this way we prove 
that, using the proper mode expansions is equivalent to working in the reduced phase space. 
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In section 6, we apply the machinery developed in the previous sections to the case of mixed 
boundary conditions, i.e. , we find the constraints chain, the Dirac bracket and the reduced 
phase space. The new and interesting result of this case is that, the Dirac bracket of two 
field components is obtained to be non-zero and hence, in the quantum theory these field 
components are noncommuting. The last section is devoted to the concluding remarks. 

2 Review of Dirac Procedure 

Given the Lagrangian L(q,q) (or L((f),d<f)) in a field theory), the Lagrangian equations of 



where Lj are Eulerian derivatives, W{j(q, q) = q^q^ is called the Hessian matrix, and «j = 
§qZ ~ ( a^dqi ) • ^ l^tfl = ^' the Lagrangian is called singular and in this case the number 
of equations containing accelerations are less than the number of degrees of freedom. Hence 
a number of Lagrangian constraints, J a (q,q) = 0, emerges (To obtain these constraints we 
should simply multiply both sides of (2.1) by the null eigenvector A" of W, so 7 a (g, q) = A"aj, 
||.). Then we should add the time derivatives of constraints, A f a (q, q), to the set of equations 
of motion to get new relations containing the accelerations. As a result two cases may happen 

1) the rank of equations with respect to acceleration is equal to the number of degrees of 
freedom. 

2) new constraints, acceleration free relations, emerging. 

In the first case the equations of motion can be solved completely, however, the solutions 
should obey the acceleration free equations, the constraints. In the second, the derivatives 
of new constraints and derivatives of previous constraints should be added to the equations 
of motion, and the same scenario should be repeated. 

At the end, there may remain a number of undetermined accelerations; it is shown that 
they correspond to the gauge degrees of freedom and are related to the first class Hamiltonian 
constraints. Moreover roughly speaking, there may exist some degrees of freedom which have 
no dynamics and are completely determined via the constraints. These are related to the 
the second class Hamiltonian constraints |J. 

Let us study the Hamiltonian formulation. Singularity of the Hessian matrix, J^-, implies 
the Legendre transformation, (q,q) —>■ {q,p), to have a zero Jacobian and hence, the set of 



motion are: 



Li = Wijq) + a>i = 0, 
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momenta, pf, 

dL , . 

* = W (2 ' 2) 

are not independent functions of q and q. So a number of Hamiltonian primary constraints 
turns up 

$< 0) (g,p) = 0. (2.3) 
It can be shown that [|TJ dynamics of any function in phase space is obtained by 

g^{g,H T } P . B ., (2.4) 

where weak equality, ~, is the equality on the constraint surface, and 

H T = H + A a $ a , (2.5) 

is the total Hamiltonian, A a being the Lagrange multipliers. 

Like the Lagrangian case the consistency conditions of the primary constraints should be 
investigated, i.e. the constraints should be valid under the time evolution: 

$(°) « {^\H T } P . B . « {*M,H} + X b {^\ $i 0) } « 0. (2.6) 

If the above relation dose not hold identically, then two possibilities remains i) {$^°\ }'s 
vanish weakly. In this case new Hamiltonian constraints 

$« = {<J>i°),#}, (2.7) 

turns up. 

ii) $1°''} do not vanish, yielding equations for determining A a . 

In general, depending on the rank of the matrix {$„ \ ( &i ' ) }, we may have a mixture of 
two possibilities. That is, some of the Lagrange multipliers are determined and a number of 
new constraints emerge. Here we do not bother the reader with the details. A complete and 
detailed discussion can be found in 0. 

Now the consistency conditions of $W should be verified which may result into some 
new constraints $^ 2 -*. The procedure goes on, and finally we end up with some constraint 
chains. Roughly speaking, each chain terminates if a Lagrange multiplier is determined or, 
if we get an identically satisfied relation. The latter case happens when the last constraint 
has weakly vanishing Poisson bracket with the primary constraints and the Hamiltonian. 

We denote the set of constraints $*- 2 ), .... as secondary constraints . These are really 
consequences of primary constraints while the primary constraints, by themselves have their 
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origin in the singularity of the Lagrangian (singularity of the Hessian matrix). In a pure 
Hamiltonian point of view, however, the origin of primary constraints is not of essential im- 
portance. In any way given some primary constraints, we should build the total Hamiltonian, 
(2.5), and check their consistency. 

There is another important classification of constraints: If the Poisson bracket of some 
constraint with all the constraints in the chain vanishes, it is called first class. And if the 
matrix of mutual Poisson brackets of a subset of constraints, C MN , 

C MN = (2.8) 

has the maximally rank, it is invertible, we deal with second class constraints. It is shown that 
a constraint chain terminating with an identity, is of first class and ending with determining 
Lagrange multipliers are of second class ||. To find the dynamics of a system with second 
class constraints, one may use the Dirac bracket, 

{A,B} D . B , = {A,B} P . B . - {A,$ m }p.b. (C-y^^N^jp.B.. (2.9) 

The important property of the Dirac bracket is that for an arbitrary function A and for all 
second class constraints $m, 

{$m,A} d . b . = 0. (2.10) 

It can be shown that, using the Dirac brackets instead of Poisson bracket, is equivalent to 
priory putting the second class constraints strongly equal to zero. 

For second class constraints we can always find a canonical transformation such that 
the constraints, $m, he on the first 2n coordinates (qi, q n ',Pi, ■■■,Pn) of the phase space 
and the remaining degrees of freedom, (Qi, Qn-u] Pi, ■•■> Pn-u) are unconstrained. The 
Dirac bracket in the original phase space is equal to the Poisson bracket in the space 
(Qi, Qn-u] Pi, •••) Pn-u), the reduced phase space |], 10, |TTJ. Although finding the above 



canonical transformation is not an easy task, for the case we study in this paper, bound- 
ary conditions as constraints, we show that using the suitable mode expansions, is in fact 
equivalent to going to reduced phase space. 



3 Boundary Conditions as Constraints 

Boundary conditions are acceleration free equations which in general are not related to a 
singular Lagrangian. To visualize this point, let us take a simple (1 + 1) field theory as a toy 
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model 

S = l ( dx [ t2 dt i(d t <P) 2 - (d x <P) 2 ]. (3.i; 



2 Jo Jti 

Variation of the action with respect to <p gives 

5S= tdx t 2 dtL{4>)5(f> + [ t2 dt(d x (p)5(j)\ l + f ' db(dt4>)5<j>\%, (3.2) 
Jo Jti Jti Jo 

where L(0) = d 2 (p — d 2 <j), is the Eulerian derivative. For an arbitrary 5<p, variation of the 
action vanishes if the three terms in the above equation vanish independently. The first 
term in (3.2) leads to equations of motion and the last term to the initial conditions. 
The second term which is called the surface term, results in the boundary conditions. For 
this term to vanish, there are two choices 8(j>\ boundary = 0, Dirichlet boundary conditions, 
or d x (j)\boundary = 0, Neumann boundary conditions. The boundary conditions unlike the 
equations of motion, are acceleration-free equations and should be held at all the times. In 
other words, they can be treated as Lagrangian constraints. To clarify this point we repeat 
the above argument in the discrete version: 

1 rt2 N N ~ 1 1 

S = o / #£ e W*) 2 - E - </>m) 2 , (3-3) 
A Jtl i=0 i=o e 

<f>i(t) = (f)(x,t)\ x=x . ;x n = ne, (3.4) 

and e = so that e — » (N — >• oo) reproduces the continuum theory. 
Demanding the variation of (3.3) to vanish, leads to 

edUo = h<Pi~M, (3.5) 
ea t 2 i = ^(0 i+1 -20 i + i _ 1 ), z^O, N (3.6) 

edt<f) N = \(j) N - (pN-i)- (3.7) 

Taking the continuum limit, assuming that acceleration of the end point are finite, the 
equations for 0, N give 

lim -(<f>i — 4>o) = and lim -(4>n — 4>n-i) = 0. (3.8) 
e e 

Hence in the continuum limit equations of motion for the end points give acceleration free 

equations, the Lagrangian constraints, where as (3.6) leads to L(0) = 0, which actually 

contains the acceleration term. 



2 It is worth noting that we still have the option 6<po or 54>n = which in the continuum limit translate 
into the Dirichlet boundary conditions . 
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A new feature appearing here is that, unlike the usual Lagrangian constraints, boundary 
conditions are the constraints which are not consequences of the singularity of Lagrangian, 
but a result of taking the continuum limit. 

4 The Hamiltonian Setup 

In this section, by going to Hamiltonian formulation, we apply the Dirac procedure to a field 
theory with given boundary conditions. Again, we take our simple toy model and treat the 
boundary conditions as Hamiltonian primary constraints: 

$(°) = d x <P\ x=0 . (4.1) 

Here we explicitly work out Neumann boundary condition at one end , the Neumann bound- 
ary condition at the other end and the Dirichlet cases can be worked out similarly. The total 
Hamiltonian is built by adding the constraint to the Hamiltonian by arbitrary Lagrange 
multiplier 

H T = H + \$ {0 \ (4.2) 

with 

H=\fdx n 2 + {d x 4>)\ (4.3) 
I Jo 

n = dt<j>. (4.4) 

We should remind that as discussed in sec. 2, the appearance of the constraints (4.1) is not a 
consequence of the definition of the momenta for an ordinary singular Lagrangian and hence, 
the transformation (4.4) between the velocities and momenta is well-defined and invertible 
throughout all the points, even at the boundaries. 
Now we should check the consistency condition 

$(°) = {&°\H T }p. B . = d x U\ = (4.5) 

which leads to the secondary constraint, It should be noted that to obtain (4.5), 

although the conditions are imposed at the boundaries, the fields can safely be extended 
into neighbourhood of the boundaries and we can use = / 5(x)d x (f)dx. 
We should go further: 

$(!) = H T } = + A{$«, $(°)} = 0. (4.6) 
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The second term on the right hand side, 



A{$( 1 ),$(o)} = f 5{x)5(x'){d x Il,d xl (f)}dxdx' 

(4.7) 

= — / 5(x)5(x')d x d x '5(x — x')dxdx', 

is not well-defined, and formally can be written as d^S(x — x')\ x=x i =0 . This term compared 
to the first term is infinitely large. The only way to impose the consistency condition on the 
constraints is 

A = 0, (4.8) 

and 

{$ (1) ,#} = 0. (4.9) 

There is a new feature appearing which is not any of the cases i) and ii) discussed in section 
2. The consistency condition, (4.6), reduces to two equations, (4.8) and (4.9), and although 
the Lagrange multiplier is determined the constraint chain is not terminated. 

The above discussion can be better understood if the calculation is regularized by con- 
sidering the discrete case. Using the the discrete version of equation (4.6), A is turned out to 
be of the order of e, going to the continuum limit is vanishes, and the other term, H}, 
should vanish separately. 

Defining H} as & 2 \ the other secondary constraint, we find 

$< 2 > = fl£4,. (4.io) 

We should go on: 

$ (3) _ £( 2 ) = {$V) jHt } = {$(2), #} = flgn|o. (4.11) 
This process should be continued and finally we are left with an infinite number of constraints: 



4 n+1 V|o n = 0,2,4,... 

(4.12) 

0( n >II|o n = 1,3,5,... 



Exhausted the constraint consistency conditions, we show that the Poisson bracket of the 
constraints, 

C mn = {$M, (4.13) 
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is non-singular and hence, the set of constraints (4.12) are all of second class. To show this 
first we calculate 



m,n= 0,2,4, ... 

m,n = 1,3,5,... ( 4 - 14 ) 

/ 5(x)S(x')d^ +1 d^8(x - x')dxdx' m = 0, 2, 4, n = 1, 3, 5, ... 



To find detC, the non-zero elements should be regularized. This regularization can be done 
by two methods, discretization or using a limit of a regular function, e.g. the Gaussian 
function, to represent 5(x). Here we choose the second, but one can easily show that the 
other method gives the same results. Inserting 



S(x — x) = lim 



(4.15) 



into (4.14) we find 



J 5{x)5{x')d™ +1 d™,5{x - x')dxdx' = ^J e -(m+n+2) H m+n+1 (0) 
= ^i(_2)( n + m + 1 )/ 2 e-( m + n+2 )(m + n)!!, m = 0,2,4, ...,n = 1,3,5, 



(4.16) 



H n (0) denotes the Hermite polynomials at x = |H^| . Putting these together, C is finally 
found to be 

C = A®B, (4.17) 

where 



A 



V 



-1 



and B is an infinite dimensional matrix with 



B n 



4<- 2 > 



(n+m— 1) 



1} ((2(m + n) -3)!! 



r -2(m+n)-l 



m, 72, = 1,2 



(4.18) 



It is straightforward to show that the matrix B has non-zero determinant, i.e. the matrix C 
is invertible, and hence all the constraints in the chain are second class. One way to consider 
all of them is using the Dirac bracket. To find Dirac bracket of any two arbitrary functions 
in the phase space, it is enough to calculate Dirac brackets of (0, 0) , (0, II) and (IT, IT) 



! Find more detailed calculations in the appendix. 
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{<j>{x)A{x')} D . B . = -{4>{x)M m) }C;^ (n) A{x')} = o. (4.19) 
{u(x),u(x')}d.b. = -{n(x),$( m )}c^{$w,n(x')} = o. (4.20) 

{<P(x),U^')}d.b. = 5(x-x')-{<P(x),<S>^}C m 1 n {& n \ll(x>)} 

(4.21) 

= S(x-x')-R(x,x'). 
Without using the explicit form of C^ 1 one can show 

R(x,x') = Ke5(x)5(x'), (4.22) 

where k is a numeric factor. To find k, let us obtain the Dirac bracket of the constrain $° 
with an arbitrary function /, using (2.10) we should have 

R0(x)| o ,/(0,n)} D . B . = J 5(x)d x {<f>(x),f} D . B . = 0. (4.23) 
Denoting d ^ x ,^ = g{x'), we can write 

J 8{x)d x {4>{x),Tl{x')} D . B .g{x') = 0. (4.24) 
Inserting (4.21) and (4.22) into (4.24) reduces to 

J (dj(x) + Ke5(x)dJ(x))g(x) = 0. (4.25) 
Remembering (4.15), we find 

k = -ypH. (4.26) 

Hence 

{(p(x), U(x')} d .b. = 5(x - x) + Ke5(x)5(x'). (4.27) 

Appearance of regularization parameter, e, in the Dirac bracket sounds bad, but since the 
second term has two delta functions, to be of the same order of the first term, in fact an e 
factor is necessary. We will clarify and discuss this point in the next section. 

The Dirichlet boundary condition can be worked out similarly. In this case the constraint 
chain is obtained to be 

0("V| O n = 0,2,4,... 



$0) 



(4.28) 

a^nio n = 1,3,5, ... 



Performing the calculations, one can show that the Dirac brackets are like the Neumann 
case, except for the k factor, which is 
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5 Mode expansion and Reduced Phase Space 



In the previous section we showed that a field theory subjected to the Neumann or Dirichlet 
boundary conditions is a system constrained to an infinite chain of second class constraints. 
As mentioned in sec. 2, for a system with second class constraints, there is a subspace of 
phase space which is spanned by a set of unconstrained canonical variables, the reduced 
phase space. The important property of these variables is that, Poisson bracket in terms of 
them is equivalent to the Dirac bracket defined on the whole constrained phase space. 

In this section we will explicitly find the reduced phase space and show that it is in fact 
equivalent to phase space determined by the Fourier modes. 

Let us consider the Fourier transformed variables 



(5.1) 



<f>(x) = ^7= j (j)(k)e ikx dk , (f)(k) = -j= J <j)(x)e~ ikx dx 
n(x) = ^ J U(k)e~ ikx dk , II(fc) = ^ J U(x)e ikx dx. 
One can easily show that the above transformation is canonical: 

{<j)(k)A(k')} = o , {n(fc),n(A/)} = o , {<j>(k)Mk')} = 5(k-k'). (5.2) 

The Neumann (Dirichlet) constraint chain, (4.12) and (4.28), in terms of the new variables 
are easily obtained: All the odd (even) moments of (f>(k) and H(k) are zero. The most general 
solution to these conditions is that (f>(k) and II(fc) are even (odd) functions of k. Then (5.1) 
gives [] 

(j)(x) = —j= I 4>{k) cos kxdk , Yi{x) = —= / H(k) cos kxdk. (5.3) 
yTT J V n J 

The main advantage of the Fourier modes, (j){k) and U.(k), is that although they are limited 
to even (odd) functions, are still canonical variables; in contrast with the original fields <fi(x) 
and U(x) which lose their usual canonical structure due to constraints. 

To compare the Dirac bracket results with those of reduced phase space, we work out 
Poisson brackets of and II (a;). Using (5.2) and (5.3) we have 

{0(x),</>(x')} = O, 

{n(x),n(x')} = o, (5-4) 

{4>(x),U(x')} — - J cos kx cos kx'dk = 5n(x,x'), 



'For Dirichlet case Cosine should be replaced be Sine. 



11 



for Neumann boundary conditions. For Dirichlet case only {0,11} differs from above: 

{<p(x), n(x')} = — / sin kx sin kx'dk = <Wx, x). (5.5) 

7T J 

Performing the integrations we have 



Sn(x, x') = S(x — x') + 8(x + x'), 
Sd(x, x') = 5(x — x') — 5(x + x'). 



(5.6) 



If we consider only the positive x's, x > 0, Sjy and 5d for x, x' 7^ are exactly 8{x — x'). For 
x, x' = 0, we regularize delta functions by 



5{x — x') + S(x + x') ■ 
5{x — x') — 5{x + x') 



2 

= 0. 



clt X — X — 



(5.7) 



Hence 5n and 8d for x > is in exact agreement with the Dirac bracket results obtained 
in previous section. The above argument clarifies, why using the usual mode expansions 
to quantize a system with Neumann or Dirichlet boundary condition, i.e. , imposing the 
boundary conditions and then quantizing, works. 



6 Mixed Boundary Conditions, Another Example 

In this section we handle a more general family of boundary conditions, mixed boundary 
conditions, which are combinations of Neumann and Dirichlet cases. It has been shown that 
these boundary conditions lead to unusual results in the context of string theory H, []]. |5[ 

As a toy model for a field theory resulting in the mixed boundary conditions consider 

S = \fdx f 2 dt[{d t ct>{) 2 - (d x <p t ) 2 + I^iho.O,,,;. (6.1) 
Z J J 1 1 

where i,j = 1,2 and is a constant antisymmetric background. Varying S with respect 
to <pi, gives: 

d 2 t 4>i - d 2 x <j>* = 0, (6.2) 
d x (pi + Fijdt^j = at x = 0, 1. (6.3) 
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Equations (6.3), as discussed in section 3, give the Lagrangian constraints. In the discretized 
version, (6.3) are the equations of motion for the end points and in the continuum limit, 
the acceleration term disappears. It is worth noting that (6.3) reproduce the Neumann and 
Dirichlet boundary conditions for F = and oo respectively. 

Now to apply the Dirac method, we go to Hamiltonian formulation: 

Ilj = d t (pi + Fijd x <f>j, (6.4) 

H=^j\u i - Fijd^j) 2 + (d x <p t ) 2 dx, (6.5) 
and the primary constraints, 

= Mx)\x=*>, (6-6) 

with 

$i(x) = Mijd x (f)j + Fijllj = , Ma = (1 - F 2 )ij. (6.7) 

Note that in this case the Lagrangian constraints, (6.3), depends on velocities, and as men- 
tioned before, the transformation (6.4), is non-singular and the Lagrangian constraints can 
be translated into Hamiltonian constraints, (6.7), without any difficulty The consistency of 
the primary constraints should be verified: 

= {*?\H T } = {$<°U} + A,{$f\$f } = 0. (6.8) 

The first term is easy to work out: 

= {$f U} = W =0 . (6.9) 

Similar to the arguments of sec. 4, {<&j°\ is infinitely large compared to the first term, 
and the only way for (6.8) to be satisfied is 

Xi = and = 0. (6.10) 



Again, although the Lagrange multiplier, Aj, is determined, there are secondary constraints, 
$^ = 0. Moreover we have the advantage that Aj disappears in the remaining steps. 
Direct calculations on the consistency conditions for the constraints leads to the chain 

d^i\o n = 0,2,4,... 



(6.11) 

4")ILJ n =1,3,5, ... 
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To verify that these constraints are really second class, we study the matrix, C™ 
{$("0, $(»>}: 



mn 



m, n = 1, 3, 5, ... 

-2(MF) ij / 5(x)5(x')a™ +1 <9^5(a; - x')dxdx' m,n = 0, 2, 4, ... 

My / - x')dxdx' m = 0, 2, 4, n = 1, 3, 5, ... 

C can be written in the form of 



(6.12) 



where F is a 4 x 4 matrix: 



C = F®B, 



^ -2{MF) M ^ 



V 



-M 







(6.13) 



(6.14) 



/ 



and B given by (4.18). In section 4, we discussed that B is invertible. Since det F ^ 0, C 
is invertible too, hence all the constraints in the chain (6.11) are second class. 
One can show that the fundamental Dirac brackets £1X6 cLS following 



(6.15) 

(6.16) 
(6.17) 



= (-2M' 1 F) lJ (e 2 ^5(x)5(x')) 

{u t (x),u,(x')} D . B . = -{Mx^KC- 1 )^^,^')} = o. 
{j> i (x),Yi j (x , )} D . B . = s(x-x') - {Mx)^ ) }(c- 1 )r{®i n \n j (x')} 

= 8(x — x') — R(x, x 1 ) = 6n(x, x'). 

The important result of the mixed case is (6.15); the Dirac bracket of two field com- 
ponents are non-zero. This means that in the quantized theory these field components are 
noncommuting. In the string theory, where the fields describe the space coordinates, (6.15) 
tells us that, the space probed by open strings with mixed boundary conditions is a non- 
commutative space 0, H, |] . 

Using the canonical (or Fourier) transformations, (5.1) and (5.2), we can explicitly build 
up the reduced phase space for the mixed case. Let $i(k) represent the Fourier modes of 
defined in (6.7), 



y In J v 2n J 



(6.18) 
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Using (5.2), Poisson brackets of $j(fc) and Hi(k) can be worked out. Imposing the constraints 
(6.11), we find that &i(k) and Hj(k), are odd and even functions of k respectively: 

$j(x) = -^= [ §i(k)smkxdk , UAx) = -^= / UAk) cos kxdk. (6.19) 

a/7T J \/ 7T J 



Remembering (6.7), we can derive the field components: 

Mr. 1 r -dk 



<j>i{x) = —*Lr / —— (^j(k) cos kx + F jk U k (k) sin kx), (6.20) 

\/7T J k 



which explicitly satisfies the mixed boundary conditions. 

Derived the mode expansions of the fields and their conjugate momenta, we finds their 
Poisson brackets: 

= \S fl£[( M ~ lF )*{Mk), Mk')}M^ cos kx sin k'x'+ 
+ (M- 1 F) jk {Tl k (k),<S> l (k')}M u 1 cos k'x' sin kx] 
= ~V / ^-(M^ 1 F) ij (cos kx' sin kx + cos kx sin kx') (6-21) 
= (M~ l F)ij J x (S N (y, x') - 5 D (y,x'))dy 
= -2{M- l F) l3 j x 5{y + x')dy. 

Since for x, x' > 

{2 x — x ' — Q 
(6.22) 
otherwise, 

(6.21) is non-zero only for x, x' = 0: 

{<f> i (0),<f> j (0)} = -2(M- 1 F) ij . (6.23) 

Comparing (6.21) and (6.15), we find that they are exactly the same. In other words, (6.19) 
and (6.20) are functions defining the reduced phase space. 

In the context of string theory, (6.21) implies that the end points of open strings subjected 
to mixed boundary conditions are living in a noncommutative space. The mixed open strings 
appear when we are studying D-branes in a NSNS two-form background. In this case, (6.21) 
tells us that the world-volume of such branes are noncommutative planes. 

We can also calculate {IL(x), IL(x')} and {<pi(x), Uj(x')}. The results are in exact agree- 
ment with (6.16) and (6.17). 
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7 Concluding Remarks 



In this paper, we studied the old and well-known problem of field theories with boundary 
conditions from a new point of view. We discussed that in the Lagrangian formulation 
boundary conditions are Lagrangian constraints which are not a consequence of a singular 
Lagrangian. For further study we built the Hamiltonian formulation, and considered bound- 
ary conditions as primary constraints. Asking for the constraints consistency conditions we 
found two new features in the context of constrained systems 

1) Although the Lagrange multiplier in the total Hamiltonian is determined, the con- 
straints chain is continued. 

2) Boundary conditions are equivalent to an infinite chain of second class constraints. 
Constructing the Dirac brackets of the fields and their conjugate momenta for these 

second class constraints, we showed that the method based on mode expansion, is equivalent 
to working in the reduced phase space. 

The relation between Hamiltonian method we built here and the usual method of im- 
posing boundary conditions in the equations of motion, can simply be understood. In the 
former, to ensure that boundary conditions are satisfied, we make the Taylor expansion of 
boundary conditions as a function of time, and put all the coefficient equal to zero. These 
coefficients are exactly our constraint chain. But, in the latter, the Fourier mode expansion 
is used and boundary conditions are guaranteed by choosing all the Fourier components to 
satisfy boundary conditions. 

In the last section of the paper, we handled the mixed boundary conditions which seems 
to be an exciting problem in the context of string theory [0]. Having noncommuting field 
components, is the interesting feature appearing in this case. Besides the string theory, 
mixed boundary conditions can be encountered in the context of electrodynamics having an 
extra #-term: 



In the above action 9 plays a role similar F in our toy model. Varying the action gives 
a surface term, vanishing of which leads the mixed boundary conditions. Quantizing this 
theory is an interesting problem we postpone it to future works. 
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Appendix: In this appendix we present some of the calculation details 



L^ m \x) = {n(a;),$ (m) } 



m = 0,2,4,... 

{(f>(x), a^n( m )l°} = k m (x) m = 1, 3, 5, .... 

{n(x),9™^ m )l°} = l m (x) m = 0,2,4,... 
m = 1, 3, 5, .... 



k m {x) = [ d%8(x - x')S(x')dx' = -L e xp(^) -^H m (0) = 5(x)k m . 

r( x ) = - [ d™ +1 6(x' - x)8{x')dx' = -^exp(^) -±-H m+1 (0) = 5(x)k m+1 . 
where, H m (0) is the Hermite polynomial at zero. Then one can easily work out 

{<t>(x),n{x')} D . B . 

{4>(x),II(x')} d .b. = S(x - x) + k m+ ik n B~ 1 n 5(x)S(x') 

The power of e in km+iknB^, can be read off from the explicit form of k m and B mn , and 
the results is k m+ ik n B^ n = ke. 

Calculations for the mixed boundary conditions can be performed similarly. 
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